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FOREWORD 


This  report  is  submitted  to  the  Ballistics  Research  Laboratory  in  partial  fulfillment  of 
Delivery  Order  0001  of  Contract  No:  DAAA15-90-D-1Q02. 

The  BRL  Project  Officer  is  Mr.  Richard  Pearson.  The  SPARTA  Program  Manager  is  Mr. 
Gregory  Mason.  Mr.  Daniel  Nowlan  performed  the  elegant  control  system  dynamic 
analysis  and  Dr.  Irving  Osofsky  provided  valuable  technical  advice. 
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PREFACE 


On  30  September  1992,  the  U.S.  Army  Ballistic  Research  Laboratory  (BRL)  was 
deactivated  and  subsequently  became  part  of  the  U.S.  Army  Research  Laboratory  (ARL)  on 
1  October  1992. 
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1 .  INTRODUCTION 


The  US  Army  Ba"istic  Research  Laboratory  (BRL)  is  currently  modifying  an  existing  shock 
tube  located  •=>!  the  Aberdeen  Proving  Grounds  to  demonstrate  Large  Blast/The rma! 
Simulate-  fechno!ogies  and  to  provide  high  fidelity  air  blast  environments  for  nuclear 
effects  testing.  This  facility  will  be  the  first  large  shock  tube  to  use  heated  driver  gas  to 
acmeve  desired  air  blast  waveforms  m  the  test  section.  The  driver  gas  must  be  heated 
and  the  driver  quickly  pressurized  to  minimize  heat  loss  to  the  driver  walls;  this  required 
innovative  solutions  to  pumping  hot  gas. 

Under  BRL  sponsorship  SPARTA  developed  and  demonstrated  a  driver  filling  method 
which  pumps  liquid  nitrogen  (LN)  through  a  previously  heated  Pebble  Bed  Heater  (PBH) 
thereby  vaporizing  the  liquid  and  raising  its  temperature  to  the  desired  value  in  one  pass 
(Figure  f).  A  bypass  system  allows  precise  control  of  the  output  gas  temperature  by 
selectively  mixing  LN  with  the  heated  gas  exiting  the  Pebble  Bed  Heater.  This  approach 
has  the  advantages  that;  pumping  a  liquid  is  more  efficient  than  pumping  a  gas  (if  indeed 
pumping  a  hot  gas  can  be  done  at  all),  the  LN  pump  is  much  smaller  and  much  more 
robust  than  a  gas  compressor  system  and  the  constant  displacement  pump  mass  flow 
rate  is  independent  of  back  pressure. 

SPARTA  installed  a  22  ton  Pebble  Bed  Heater  working  unit  at  BRL  and  successfully 
demonstrated  its  performance  (Reference  1).  Manually  operated  valves  were  used  to 
route  the  liquid  nitrogen  to  the  Pebble  Bed  Heater  in  these  initial  tests.  However,  an 
automatic  control  system  is  preferable  to  manual  operation  for  safety,  precise  gas 
temperature  control  and  efficiency  of  operation. 

1.1  Objectives 

The  objectives  of  the  present  study  were  to  design  an  automatic  control  system  which 
meets  established  requirements,  analyze  its  performance,  prepare  drawings  and  provide 
a  system  cost  estimate. 

1 .2  Requirements 

Four  driver  gas  design  conditions  were  established  by  BRL  (Table  1).  The  driver  filling 
strategy  is  to  pump  a  constant  temperature  gas  (at  or  above  the  design  temperature)  until 
the  driver  gas  reaches  the  design  pressure  (Reference  2).  During  the  filling  process,  the 
PBH  back  pressure  will  rise  from  ambient  to  the  peak  value  approximately  linearly 
(depending  on  the  magnitude  of  the  heat  loss  to  the  walls). 

1 .3  Scope 

Valves  and  valve  properties  were  selected  from  vendor  supplied  information.  Appropriate 
valve  settings  were  established  based  on  a  thermal  hydraulic  model  which  considered 
pressure  drops  in  the  primary  and  bypass  paths.  Control  system  response  was 
established  based  on  analytical  models  of  the  valve/actuator  motion.  Cost  estimates  were 
based  on  vendor  quotes  and  engineering  estimates  by  experienced  personnel. 
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2.  GAS  SUPPLY  CONTROL  SYSTEM  DESCRIPTION 


The  control  system  analysis  starts  with  a  control  system  layout  and  performance 
characteristics  of  the  specific  valves  used  to  direct  the  flow  to  the  PBH  and  the  bypass 
system. 

2.1  Control  System  Layout 

Liquid  nitrogen  is  pumped  from  a  tank  through  a  series  of  pipes  and  valves  to  the  Pebble 
Bed  Heater.  Numerous  diagnostic  measurements,  exhaust  valves,  check  valves  and 
control  valves  are  used  to  control  the  fluid  flow  as  indicated  in  Figure  2  which  is  taken 
from  Reference  3.  The  PBH  control  valves  denoted  as  V7  and  V8  are  of  specific  interest 
to  this  study. 

The  PBH  controls  perform  the  functions  of  regulating  the  outlet  temperature  of  the  PBH 
mixer  to  a  specified  setpoint.  The  PBH  operates  by  receiving  LN  from  the  high  pressure 
pumping  system  and  branching  the  LN  flow  into  two  subsystems,  the  pnmary  pebble-bed 
and  the  bypass  thermal  mixer.  Precise  control  of  valves  7  and  8  in  the  primary  and 
bypass  lines  is  required  to  produce  a  stable  outlet  temperature  of  the  PBH. 

2.2  Control  Valves 

Globe  valves  have  been  selected  for  the  PBH  valves  because  of  their  rugged 
construction,  cryogenic  rating  and  high  capacity.  Valve  position  control  is  achieved  by  the 
use  of  a  closed  feedback  contro!  loop  to  the  actuators  using  the  outlet  temperature  of  the 
mixer  as  the  feedback  sei  ,sor  point.  Valtek  Mark  One  valves  and  Linear  Spnng  actuators 
were  chosen  to  develop  valve  performance  characteristics  and  establish  cost  estimates 
(Reference  4). 


TABLE  1  DRIVER  GAS  DESIGN  CONDITIONS 


CASE  NUMBER 

MAXIMUM 

MAXIMUM 

GAS 

GAS 

BACK 

BACK 

TEMPERATURE 

TEMPERATURE 

PRESSURE 

PRESSURE 

(MPa) 

(ATM) 

(K) 

(R) 

1 

12.8 

129 

663 

1193 

2 

7.8 

79 

468 

842 

3 

2.9 

29 

361 

650 

4 

0.98 

10 

288 

.. 

518 

LNBP 


Figure  2.  Gas  supply  system  schematic. 


3.  THERMAL  HYDRAULIC  ANALYSIS 

Engineering  procedures  were  developed  from  analytical  models  based  on  conservation 
of  fluid  momentum  and  energy.  All  Pebble  Bed  Heater  thermal  hydraulic  processes  are 
relatively  slow  so  transients  are  not  important. 

3.1  Model 

The  relationship  defining  the  temperature  of  the  nitrogen  gas  leaving  the  thermal  mixer 
is 


rt>T(Cpnhv)  ~  rh-, ( Cp  Tpg  +hv)  + rh8 Cp TLN 


EQ(1) 


Conservation  of  mass  gives 


mT  -  m7  +■  ma  EQ  (2) 

Substituting  (2)  into  (1)  and  solving  for  temperature  gives 

(  Tpb-Tl^£.)  EQ  (3) 

mT  Cp 

For  convenience  we  define  a  reference  temperature 

V  T„-  rL„+  £  EQ  (4) 

Both  rti7  and  m8  are  a  function  of  time  during  the  operation  of  the  PBH  due  to  changing 
back  pressure  and,  near  the  end  of  the  run,  changing  PBH  exit  temperature.  Control 
valves  located  on  the  primary  PBH  supply  line  (V7)  and  the  mixer  bypass  supply  line  (V8) 
are  assumed  identical  with  respect  to  flow  capability. 

The  mass  flow  rates  are  determined  by  equating  the  pressure  drops  in  the  primary  and 
bypass  legs.  In  the  primary  system  the  pressure  drops  are  due  to  the  pipes,  valve  7,  the 
elbows  and  the  PBH;  in  the  bypass  system  the  pressure  drops  are  due  to  the  pipes  and 
valve  8.  Following  Reference  5  we  have  for  the  primary  leg 

Pipes 


DELP1  =  P  -L  fpl l  =  Cl  U72  EQ  (5) 

2  p  D  7 
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Elbows 


DELP2  -  N  JL^l  fM  =  C2  U7  2  EQ  (6) 

Valve 

.  A  ,  l2  f  >2 

DELP3-  P  U’A  4  -  C3  JL  U  >  EQ  (7) 

CV  x  4  7 

v  / 

Pebble  Bed  Heater 

DELP4  -  f„  -  C4  Up2  EQ  (8) 

Here  the  pressure  drop  constant  was  selected  such  that  the  pressure  drop  equalled  20 
psi  based  on  BRL  measurements  on  the  existing  PBH. 

For  the  bypass  leg 

Pipes 

DELP5  -  i-U±  fp  h.  -  Cl  U  2  EQ  (9) 

2  D  8 

and 

Valve 

DELP6  -  P  U*  A  ^  =  C3  (  JL  ]  U2  EQ  (10) 

CV  x  '  8 

L  J  V  ^  / 

Summing  pressure  drops 

2  /  \-2 

+  C4  ]  -  U2  [  C5  +  C6  JL  ]  EQ  (It) 

l  ^ 

Equation  1 1  can  be  rewritten  as 


U2  [Cl  +  C2  +  C3  _ 

^7 
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U72  C7  (x)  -  U2  ca  (x) 


EQ  (12) 


and  therefore 


Us 

U7 


(  \ 

C7 

ca  j 


1/2 


EQ  (13) 


where  it  is  understood  that  the  terms  C7  and  C8  are  functions  of  x. 

Since  the  liquid  nitrogen  is  essentially  incompressible,  Equation  3  can  be  rewritten  as 


T-  Tpg- 


U 


(U7+  U6) 


Tpg  (  1 


f  V/2 

C8 

C7 


)"  tr 


EQ(14) 


which  gives  the  output  gas  temperature  as  a  function  of  valve  position.  Setting  the  output 
temperature  at  the  desired  or  set  point  control  temperature 

T  =  Tr 


and  performing  a  bit  of  algebra,  there  results 


(  \ 

x 

\  P 


C3 


(  T  > 

2 

7 - ^ 

-z 

T«  -1 

(  C5  +  C6 

X 

) 

T  -  T 

{  IPB  >C  j 

10 

)  -  (  Cl  +  C2  +  C4  ) 


V2 


EQ  (15) 


3.2  Calculations 

Sample  calculations  are  made  to  indicate  nominal  valve  settings  and  their  sensitivities. 
The  basic  procedure  is  to  set  one  valve  (e.g.,  V7)  at  a  single  setting  and  control  the  PBH 
output  temperature  with  the  other  valve  (V8).  It  is  desirable  that  the  valves  be  operated 
in  mid  range,  if  possible,  simply  to  avoid  fine  settings  and/or  slamming  into  the  stops 
unnecessarily.  The  design  test  conditions  are  addressed  first  at  their  peak  back 
pressures.  Then  the  effect  of  back  pressure  is  considered. 

3.2.1  Design  Test  Conditions 

A  broad  range  of  PBH  output  temperatures  are  achievable  if  the  PBH  is  heated  to  2000 
°R  (1110  °K).  The  peak  design  temperature  of  1193  °R  (663  °K)  is  achieved  with  a 
secondary  valve  relative  opening  of  0.2  to  0.5  as  the  primary  valve  relative  opening 
ranges  from  0.25  to  1.0  (Figure  3). 
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BRL  PBH  OPERATION 

MIXER  OUTLET  TEMPERATURE 


8  FRACTION  OPEN 


Considering  each  design  case  individually,  the  sensitivity  to  outlet  temperature  is 
established  by  calculating  the  required  valve  settings  for  the  nominal  design  condition  and 
for  the  next  lowest  outlet  design  condition  (Figures  4  to  7).  Initial  bed  temperatures  were 
picked  based  on  enthalpy  scaling  from  present  test  results  (i.e.,  thermal  energy  required 
is  proportional  to  the  mass  and  temperature  of  gas  required).  However,  while  technically 
feasible  to  operate  at  PBH  temperatures  near  the  desired  outlet  temperature  for  the  lower 
pressures,  a  minimum  bed  overheat  of  500  °R  was  evolved  by  trial  and  error  to  insure 
robust  valve  control  (see  for  example,  Figure  8  compared  to  Figure  6  and  Figure  9 
compared  to  Figure  7). 

3.2.2  Effect  of  Back  Pressure 

All  of  the  pressure  drop  mechanisms  considered  except  the  PBH  assumed  that  the 
nitrogen  was  liquid;  therefore,  only  the  PBH  pressure  drop  will  be  a  function  of  back 
pressure.  We  have 


DELP4  -  f  {Re)  £  UPg  EQ  (16) 


The  Reynolds  number  range  considered  is  about  1 0  to  1 000  and  according  to  Reference 
5,  the  pressure  drop  across  a  porous  media  is  not  a  strong  function  of  Reynolds  number 
in  this  range.  Further,  the  mass  flow  is  constant  and  the  nitrogen  gas  phase  follows  the 
perfect  gas  law  which  results  in 


DELP4  -  I 
P 


EQ  (17) 


Thus,  at  a  given  bed  temperature,  the  pressure  drop  is  inversely  proportional  to  the  back 
pressure.  Figures  1 0  and  1 1  indicate  valve  setting  combinations  required  at  a  lower  back 
pressure  for  each  of  the  elevated  design  temperature  conditions.  While  there  is  some 
effect  (e.g.,  the  required  V8  relative  opening  is  lowered  about  0.1  for  a  70  percent  V7 
opening),  the  required  valve  positions  are  well  within  the  operating  capability  of  the  gas 
supply  system  and  the  long  pumping  cycle  (order  of  minutes)  allows  plenty  of  time  for  the 
automatic  control  system  to  adjust. 
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BRL  PBH  OPERATION 

REQUIRED  VALVE  OPENING 
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■ening.  TPBH  =  2000  R.  PBACK  =  130  ATM. 


BRL  PBH  OPERATION 

REQUIRED  VALVE  OPENING 
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Figure  5.  BRL  PBH  operation  required  valve  opening.  TPBH  =  1400  R,  PBACK  =  80  ATM. 


BRL  PBH  OPERATION 

REQUIRED  VALVE  OPENING 


Fi9ure  6-  BRL  PBH  operation  required  valve  opening.  TPBH  =  850  R,  PBACK  =  29  ATM 
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R9ure  7 •  BRL  PBH  operation  required  valve  opening,  TPBH  =  600  R.  PBACK  =  10  ATM 


BRL 

REQ 


BRL  PBH  OPERATION 


BRL  PBH  OPERATION 

REQUIRED  VALVE  OPENING 


16 


0.8  T 


00  CM 

o  o 


N3d0  NOIlOVttd  L  3A1VA 


17 


0.1 -F 


4.  CONTROL  SYSTEM  DYNAMIC  ANALYSIS 


A  dynamic  model  of  the  valve  operation  and  control  system  was  developed  to  assess 
automatic  valve  performance. 

4.1  Model 

Once  the  primary  valve  has  been  set,  the  mass  flow  rate  through  the  bypass  valve  is 

ma-MFR^l l  EQ  (18) 


where  MFR  =  mass  flowrate  with  the  valve  fully  open  {back  pressure  dependent). 

The  valve  seat  is  positioned  by  a  spring  loaded  pneumatically  actuated  piston.  Modern 
globe  valves  are  fast  acting  and  highly  damped  so  the  actuator  is  modeled  as  a  damped 
spring  mass  system.  The  force  balance  equation  for  the  pneumatic  actuator  is 

mx+  /cx-  F{t)  EQ{19) 

Defining  a  valve  position  forcing  function  as 


X,  (t) 


m 

k 


and  defining  the  following  constants 


con  - 


;  - 


f 

2-i  km 


Equation  1 9  becomes 


x+2C,o)nx+  u>2n  *  -  X,(t)  Mn 


EQ  (20) 


EQ  (21) 


EQ  (22) 


EQ  (23) 


The  Laplace  transform  of  EQ  23  is 
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x(s)  - 


X,(s)  o>2n  +  (s^frQxfl) 

S2+  2  £  Cl 3n  S+  (On 


f 

X,(s)  + 

\ 


(s+2Ct aJx(O)] 


G(s) 


EQ  (24) 


where  X,(s)  is  the  transformed  valve  position  forcing  function,  x(0)  is  the  valve  actuator 
initial  position  determined  from  Figures  4  to  7  and 


G(s) 


(On 


S2+  2£(0„S+ (On 


EQ  (25) 


Since  pumping  time  (minutes)  is  long  compared  to  the  valve  response  time  (  ~  seconds) 
and  valve  damping  is  high,  a  proportional  control  system  was  selected.  Setting  the 
actuator  forcing  function  proportional  to  the  difference  between  the  measured  output 
temperature  and  the  control  or  setpoint  temperature 

X,  -  K  (T-  Tc)  EQ  (26) 


Substituting  (EQ  18)  into  (EQ  3)  gives 


T-rPS-!^-  t„  *(o 


mTL 


EQ  (27) 


Taking  the  Laplace  transforms  of  (EQ  26)  and  (EQ  27)  and  combining  with  (EQ  24)  gives 


T(s)  -  Tpg(s)  ~ 


fMFR  T  ' 
mT  L 

v  r  ) 


G(s) 


K(T(s)~  T As)) 


(s-2Ca>„)x(0) 


EQ  (28) 


or 


T  (s)  -  Tpg(s)  +  K ,  G(s) 


TASS)  -  T(s)  - 


(s+  2C<on)x(0) 


K(o2 


EQ  (29) 


where  the  nondimensional 


K^K  MfR.  Tf 

1  mrL  ' 


EQ  (30) 


Solving  for  T(s)  in  (EQ  29)  gives 
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T(s) 


Tpg  (s)  *  K,  G(s)  Tc(s)  K,  (s  +  2£wn)  x(0)  G(s) 


EQ  (31) 


UK,G(s)  Kain2(1  +  K,G(s)) 

Because  Tc  (the  setpoint  temperature)  is  a  constant  value,  the  Laplace  transform  is 


Tc  (s)  -  -2 
s 


EQ  (32) 


The  temperature  of  the  outlet  gas  of  the  PBH  (primary  flow),  Tpg,  entering  the  thermal 

mixer  can  also  be  considered  a  constant  because 

1 .  This  temperature  will  be  maintained  at  a  nearly  constant  value  until  the  pebble-bed 
temperature  at  the  end  (bottom)  of  the  pebble-bed  starts  dropping.  This  typically 
would  occur  when  approximately  80  %  of  the  process  flow  has  been  used. 

2.  During  the  last  20%  of  the  flow  period,  the  gas  temperature  at  the  outlet  will 
decrease  slowly  compared  to  the  response  rate  of  the  control  valves. 

Thus,  the  Laplace  transform  of  TPH  is 


Trsis)  - 


EQ  (33) 


Substituting  (EQ  32)  and  (EQ  33)  into  (EQ  31)  gives 


Us)  -  1 
s 


TPb*  K.  TcG(s) 


i  +  K \  G(s) 


K,(s+2&n)x{0)G(s) 
Kton2(f  +  K;G(s)) 


EQ  (34) 


and  taking  the  inverse  transform  of  (EQ  34)  gives 


T(t ) 


tpb+kjc 
1 + 


Ki(Tn-Te)  K, 


1  +  K, 


K 


x(0) 


e 


cos(co„fJ  1  +  K,  -  C2  )*  t  I— ~  —  ■  sin  (c 1  +  K:  -  i;2  ) 


EQ  (35) 


High  gain  (K,  »  1 )  is  required  to  result  in  the  steady  state  value  of  T  approaching  Tc. 
Inspection  of  EQ  35  indicates  that  for  TPaof  the  order  of  2000°R  and  Tc  of  the  order  of 
1000  °R,  K,  must  be  of  the  order  of  50  to  achieve  controlled  mixed  gas  temperatures 
within  5  percent  of  the  desired  temperature  Tc. 
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The  resuiting  solution  for  valve  position  (from  EQ  27)  becomes 


*00 

L 


rti  r  ( Tpg  7~c) 


f 

y 

1- 

\-Kzx^0) 

e"At  ( 

\ 

MFR  Tr  UK , 

U - 

j i^-C 


£ _ sin  (ajJ  UK,-?  ) 


EQ  (36) 


where  the  nondimensional 


L  1  + 

*  T*-  Tc 


4.2  Calculations 

Example  calculations  were  made  to  demonstrate  valve  performance  using  the  following 
characteristics  from  Valtek  literature: 

an  =  spring  cylinder  actuator  natural  frequency  -  2.5  cps 
£n  *  damping  ratio  *=  0.7 

Selecting  an  initial  position  of  the  primary  valve  as  70  percent  open,  the  peak  design 
condition  of  2000  °R  and  1 29  atm  back  pressure,  an  initial  valve  position  of  half  open  and 
a  K  of  50,  the  setpoint  temperature  is  achieved  in  less  than  a  quarter  of  a  second  (Figure 
13)!  Note  that  the  valve  motion  is  minimal  due  to  the  excellent  choice  of  initial  position; 
Figure  14  provides  a  better  view  of  the  motion  by  changing  the  scale  of  the  ordinate. 
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PEBBLE  BED  HEATER 

VALVE  RESPONSE 


Pebble  bed  heater  valve  response. 


PEBBLE  BED  HEATER 
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5.  Cost  Estimate 


Components  of  the  control  system  are  itemized  in  Table  2  along  with  a  suggested  vender 
and  catalog  prices.  Unburdened  hardware  costs  total  $56,000.  The  estimated  price  to 
procure,  assemble,  program,  install  and  checkout  the  system  ;s  $200,000.  Including 
hardware  purchases,  the  period  of  performance  is  expected  to  be  six  months. 
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Table  2  Automatic  Control  System  Components 
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Table  2  (Cont’d) 
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Table  2  (Cont’d) 


6.  Summary 

Conclusions  and  recommendations  reached  as  result  of  this  design  study  are  summarized 
below. 

6.1  Conclusions 

*  An  automated  control  system  for  the  BRL  t/6th  scale  shock  tube  is 
feasible  and  practical  using  off  the  shelf  hardware. 

*  Globe  valves  of  the  type  selected  for  valves  7  and  8  give  positive  control  over 
the  output  gas  temperature  with  reasonable  valve  actuator  positions. 

*  Driver  gas  backpressure  has  a  small  effect  on  control  valve  position  based  on 
the  Pebble  Bed  Heater  pressure  drop  model  developed  here. 

*  The  analytical  models  should  be  calibrated  with  test  data  for  each  design 
condition. 

6.2  Recommendations 

*  Assembly,  programming,  installation  and  checkout  of  an  automated  control 
system  should  begin  immediately. 

*  The  time  dependent  driver  gas  filling  model  developed  in  Reference  2  should  be 
coupled  to  the  gas  supply  system  quasi-steady  and  dynamic  control  models  in  a 
system  simulation  to  develop  detailed  control  strategies  for  each  test  condition. 
This  simulation  should  be  used  to  establish  test  procedures,  train  operators  and 
analyze  test  data. 

*  Pebble  Bed  Heater  and  mixer  nozzle  pressure  drop  data  should  be  obtained  for 
a  variety  of  flow  conditions  and  a  more  detailed  model  developed  for  use  in  the 
control  system  model. 

*  A  control  system  analysis  should  be  performed  for  the  existing  valves  which  are 
being  retrofitted  for  automatic  control  and  the  models  should  be  calibrated  with  test 
data. 
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Nomenclature 


A  area 

CV  valve  characteristic 
Cp  specific  heat  at  constant  pressure 
d  diameter 

f  fncticn  factor 

F  force 

G  transfer  function 

l\  heat  of  vaporization 
L  actuator  full  travel 

k  spnng  constant 

K  gain  constant 

m  mass 

N  number  of  units 

Re  Reynolds  number 

s  transform  vanable 

t  time 

T  temperature 

U  fluid  velocity 

x  position 

Greek 


to,  natural  frequency 

;  damping 

p  density 

Suberscnpts 

•  denvative  wrt  time 

second  denvative  wrt  time 
Subscnots 

7  valve  7  or  pnmary  line 

8  valve  8  or  bypass  line 

PB  pebble  bed 

p  pipe 

T  total 

90  elbow 

LN  liquid  nitrogen 

R  reference 

C  control 
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